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Abstract 

The fastest tunneling response in double barrier resonant structures is in- 
vestigated by considering explicit analytic solutions of the time dependent 
Schrodinger equation. For cutoff initial plane waves, we find that the earliest 
tunneling events consist on the emission of a series of propagating pulses of the 
probability density governed by the buildup oscillations in the quantum well. 
We show that the fastest tunneling response comes from the contribution of 
incident carriers at energies different from resonance, and that its relevant 
time scale is given by r r = nh/ \E — e|, where e is the resonance energy and 

E is the incidence energy. 
PACS: 73.20.Dx,73.40.Gk 
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The transport mechanism in resonant tunneling diodes (RTD) and the relevant time 
scale for its fastest tunneling response have been the subject of intense investigationfli'ii 
In particular, the results of the high-frequency tunneling experiments! demonstrated the 
possibility that the charge transport may occur at time scales even shorter than the lifetime 
of the quasibound state of the system. Most of the theoretical efforts to estimate the relevant 
time scales for the tunneling process have been based on stationary approaches;!@ however, it 
has been widely recognized!'!! that the analytic solution of the time-dependent Schrodinger 
equation (TDSE) provides the most reliable way to tackle this fundamental issue. Although 
there is in the literature a consensus that the buildup dynamics determines the ultimate 
speed of RTD's, the way in which the buildup governs the emission of transmitted particles 
still needs elucidation. 

In this letter we investigate the earliest tunneling events in double barrier resonant 
structures and its relation to the dynamics of the buildup process in the quantum well, 
based on the analytic solution of the TDSE. Starting from the formal solution, we show 
that the earliest evidence of transmitted carriers to the right edge of the structure (x > L) 
consists on the emission of a series of propagating pulses, which are governed by the buildup 
oscillations! in the quantum well. 

The initial condition used here is inspired on experimental situations in which the tunnel- 
ing process begins almost abruptly. It is represented by a cutoff planewave \& (x, k; t — 0) = 
0(— x)(e thx — e~ tkx ) impinging on a shutter! placed at x = 0, just at the left edge of a DB 
structure. The tunneling process begins with the instantaneous opening of the shutter at 
t = 0, enabling the incoming wave to interact with the potential, giving rise to both the 
buildup in the quantum well and the electronic transport to the right of the DB. For the 
above initial condition, the solutions ^ t (x,k;t) for the internal region (0 < x < L), and 
ty e (x, k; t) for the external region (x > L), are respectively:! 

oo 
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W = T k M{yt)-T„ k M{yt k )-i £ T n M(^J. (2) 

n=— oo 

Here 0^ = k) is the stationary wave function, T k = t(k) is the transmission amplitude, 
and the factors 4> n (x, k) = 2ku n (0)u n (x) / (k 2 — k 2 ) and T n = 4> n (L, k) exp(—ik n L) are given in 
terms of the resonant eigenfunctions u n (x) with complex eigenvalues k n = a n — ib n (a n , b n > 
0). The index n runs over the complex poles k n distributed in the third and fourth quadrants 
in the complex fc-plane. The M's are the Moshinsky functions,^ with arguments y^x^t) = 
e^ m ^{m/2ht) l / 2 [x — hqt/m], and y l q = y^(0,t), where q stands either for ±k or k± n . For the 
particular case of DB systems with isolated resonances i.e. |e n ±i — e n \ <C T n , one term is 
sufficient most of the DB diodes with typical potential parameters fall into this category. 

The emission of the earliest transmitted particles is described by the external probabil- 
ity density; this is illustrated in Fig. [I] in which we plot the normalized probability density 
|\I /6 (a;, k; t)/Tk\ 2 as a function of x for fixed values of time chosen in this example0 as t = 2 
and 10 ps. For off- resonance incidence (solid lines), the propagation of an oscillatory struc- 
ture with a sharp defined wavefront is clearly appreciated; to illustrate that this wavefront 
travels with approximately the classical speed v = (2E /m) 1 / 2 , the arrow indicates the po- 
sition x = vt, for t = 10 ps. Note that |\l/ e | oscillates around the transmission coefficient 
|Tfc| 2 , which is the expected asymptotic value as t — > oo at x = L. 

In order to analyze the tunneling mechanism at its earliest stages, we present in Fig. |2| 
a series of snapshots of this process at early times and their correlation with the periodic 
buildup oscillations observed in the quantum well. In part (a) we illustrate the buildup 
during this transient at the particular times in which the buildup in the well reaches its 
maxima and minima, as it oscillates around the stationary value (dotted line). Panel (b) 
exhibits the birth of the first pulses for the same sequence of times chosen in (a). The 
correspondence between the buildup oscillations and the successive emission of pulses is 
evident from the figure. We see that upon completion of one buildup cycle, a pulse is fired; 
in this way for example, at 2.4 ps, three cycles have been completed and up to this stage 
three pulses have already been emitted. 
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In order to find analytically how the buildup governs the fastest tunneling response and 
its relevant time scales, we shall exploit the analytical properties of the solutions. As shown 
in a recent paper,i the probability density at any position of the internal region, may be 
described by the simple formula 

¥/<f)\ 2 = 1 + e~ rnt/h - 2e- Tnt/2h cos [u n t] , (3) 

where u n = \E — e n \ /%. The oscillatory function of the above formula gives the periodicity 
of the buildup oscillations depicted in Fig. |] (a). As we shall show below this dynamical 
behavior also manifests itself outside the structure, and is the key mechanism behind the 
observed synchronization between the buildup oscillations and the pulse emission. Although 
the external solution, Eq. (|2|), is more complex than since the M functions involved in 
\l/ e depend on both position and time, a simple expression for \^/ e /T k \ 2 can be obtained. We 
start by considering the one-level formula for the external solution, 

^ = T k M(y e k ) - T_ k M(yt k ) - iT n M{ylJ - 

iT_ n M(yl_J. (4) 

From the one level expression for <p(L, A;),0 the transmission amplitude T k can be writ- 
ten as T k = 2ikexp(—ikL)u n (0)u n (L)/(k 2 — k 2 ); thus the factors iT n and iT_ n in the above 
expression can be readily identified as iT n = T k exp[i(A; — k n )L] and iT_ n = —T£ exp[— i(k — 
fc*)L]. Using the symmetry relation! for the M functions, M[y e ^) = exp(?/q 2 ) — M(— y^), 
in M(yl) and M(yf, n ), we obtain a suitable representation for \l/ e consisting on exponential 
terms and M functions with arguments of the type y e with q = —k, —k n , —A;*, which have 
vanishingly small contributions to the solution. In fact, for a fixed value of the position Xf, 
there exists a time interval starting from t > (mxf/hk), governed exclusively by the expo- 
nential terms mentioned above. From these considerations, a simple formula for \ty e /T k \ 2 
can be obtained, namely, 

| \j>e /j^j 2 = i _|_ e -r n t/h e 2b n (x f -L) _ 2 e K(x f -L) e -r n t/2h 

x cos [(a n — k)(xf — L) +oj n t] , (5) 



valid for t > (mXf/hk). By comparison with Eq. (0), we note that \^/ e /Ti t \ 2 and l^*/^ 2 
have similar time dependence. In fact, at Xf — L, reduces to 

|^ e /T fc | 2 = 1 + e- r " t/h - 2e- Tnt/2h cos [u n t] , (6) 

which coincides exactly with Eq. (Q). This analytical result confirms the qualitative dis- 
cussion of Fig. 0, that the buildup oscillations in the quantum well govern the emission of 
the propagating pulses at x = L. Therefore, both the buildup oscillations and the emitted 
pulses are characterized by the same time scales. It is straightforward to see from Eq. @ 
that the maxima of the buildup oscillations occur approximately at the time scales r m given 
by 

(2m- l)irh 



T, 



in 



AE. 



m = 1,2,3,... (7) 
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where AE n = \E — e n \ measures the deviation of the incidence energy E from resonance e n . 
Note that the same result is obtained for incidence above and below resonance, E = e n ±AE n , 
in view that Eq. (§) is irrespective to this choice. According to Eq. @, the intensity of 
the pulses reaches maximum values at x = L also at these time scales. The formation 
and emission of the first of such pulses constitutes the earliest evidence of carrier presence 
at the transmitted region; thus, the relevant time scale for the fastest tunneling response 
corresponds to m = 1, and is simply given by 

for ,„ s 

With regard to the times scales involved in the tunneling process, it is important to 
emphasize the differences between the response time r r and other relevant time scales, such 
as the lifetime, T\ = h/T, and the buildup time, T\,. In what follows we shall discuss these 
differences and their implications. 

The expressions for r r and t\ are similar since they are both inversely proportional to a 
certain energy width. However despite of this resemblance, there is a fundamental difference: 
on the one hand the lifetime is an intrinsic property of the resonant structure since it depends 



exclusively on the system parameters through the resonance width T; on the other hand, the 
time scale r r is not an intrinsic property of the system since it takes into account external 
information, namely, the energy of the incident carriers through the "off-resonance width" , 
AE. This difference is crucial for understanding the ultrafast response in DB diodes, and 
an important consequence is that depending on AE, the response time r r may be greater 
or even shorter than 77; in fact, it is easy to see from Eq. (§) that the condition for r r < 77 
is simply AE > ttT. This is illustrated in Fig. |3|, in which we compare the positions of the 
main peak of the plots of \^/ e /Tk\ 2 as a function of time and fixed position at x = L, using 
two different values of AE. 

The buildup time 77, is the duration of the transient regime in which the time- dependent 
probability density reaches its level-off given by the stationary probability density, \(f>\ . 
As shown in a recent work,! such a transient is of approximately ten lifetimes, for any DB 
structure with isolated resonances; see Fig. 2 (b) of Ref. 6. As shown in the previous 
paragraph, r r can be shorter than 77 and consequently than 77,, which means that several 
buildup oscillations (and emitted pulses) may occur before 77,. This is clearly illustrated 
in Fig. where for example at t = 2.4 ps (~ 0.38r&) significant evidence of transmitted 
carriers at the right of the structure is appreciated; in fact, up to this time the head of the 
transmitted wave has already traveled more than 10 3 nm, despite the fact that the buildup 
in the quantum well has not been fully established. 

The above results emphasize on the importance of the off-resonant carriers in ultrafast 
tunneling, which play an important role in typical tunneling experiment with RTD'sJ! where 
the energies of the available carriers are distributed within a finite interval < E < Ep, 
where Ep is the Fermi energy of the emitter. In fact, although off- resonance carriers have a 
smaller transmission coefficient, the main contribution to the tunneling current comes from 
energies near resonance rather than from the resonance alone (i.e. the Tsu-Esaki formula 
involves an integration over the whole Fermi interval, < E < Ep). As an implication of 
our analysis, when the resonance e is immersed into the Fermi sea and the resonance width 
T is such that Ep > ttT, there will be carriers of the Fermi interval fulfilling AE > ttT and 
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hence contributing with tunneling responses faster than the lifetime. The above situation 
should be present in tunneling experiments on DB systems in which Ep ^> 7rr. 

We conclude our discussion with the following remarks: (i) Since the tunneling dynamics 
critically depends on the shape of the incoming wave, whenever one deals with the problem 
of the response time, one must clearly specify to what initial condition it corresponds. In 
this respect, we stress that the response time derived here corresponds to cutoff initial 
plane waves, (ii) In view that plane waves are the "building blocks" of wavefunctions with 
more general shapes, the results obtained here may also give insight into the more intricate 
behavior expected for the case of incoming wavepackets.El (iii) The emphasis of the present 
study is the exploration of genuine quantum dynamical effects occurring in the transient 
regime; in particular, we have shown that the earliest tunneling events in a DB structure 
consist on the firing of a series of propagating pulses whose periodic emission at the right 
edge of the structure is governed by the buildup oscillations in the quantum well, (iv) The 
importance of the off-resonance incident carriers on the response time is demonstrated within 
a purely coherent picture, and a closed formula (valid near resonance) for the corresponding 
time scale has been derived from the analytic solution. 
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discussions. 
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with the first resonance. 
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FIGURES 

FIG. 1. Plot of (x,k;t) /t(k)\ 2 as a function of the distance x for two fixed values of time 
t = 2, 10 ps (solid lines). The incidence energy is below resonance E = e — AE = 74.97 meV 
(AE = ST); for comparison to the classical propagation, the arrow indicates the position x = vt 
for t = 10 ps. The special case of incidence at resonance E = 80.11 meV is also included for 
comparison (dashed lines). 

FIG. 2. Early times of the tunneling process in a DB structure. Part (a) illustrates the buildup 
oscillations inside the structure (solid line) around the stationary probability density \(f>\ 2 (dotted 
line); the potential barriers are schematically represented by the dashed-dotted lines. For the same 
sequence of times, part (b) shows the birth and emission of the first pulses which are the earliest 
tunneling events occurring in the structure. Note the synchronization with the buildup dynamics: 
for each buildup oscillation, a pulse is emitted. 

FIG. 3. Plots of (x, k; t) /t(k)\ 2 as a function of time at the fixed position x = L = 15 nm, 
for two values of AE. We see that depending on the value of AE, the response time r r may be 
greater or shorter than the lifetime of the quasibound state. 
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Figure 1 Roberto Romo. Applied Physics Letters. 




Figure 2 Roberto Romo. Applied Physics Letters. 
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Figure 3 Roberto Romo. Applied Physics Letters 




